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E ven if we use the classical physical theory of 
meteors, there may still be some possibility 
for agreement between observations and theory. 
The effect of disturbed air behind the main body 
on the motion and ablation of fragments has not 
yet been considered. This effect may have impor- 
tant consequences, observed partly as the differ- 
ence between the photometrically and the 
dynamically determined masses of the meteor 
body. By use of extreme mathematical conditions, 
this difference can be made to reach orders of 
magnitude during the latter part of the trajectory. 
However, the physical interpretation is considered 
only roughly in this paper, and the computed 
model needs further improvement. But the main 
purpose here-finding an important effect for the 
explanation of the discrepancy between the dy- 
namic and the photometric masses, especially for 
large bodies — has been achieved. 

ASSUMPTIONS 

The mathematical model used here contains the 
following simplifications: If a single meteor body 
with velocity v and radius R deposits its mass 
only in small fragments of radius R/ f the ablation 
of a spherical stony parent body is given by the 
following equations : 


dp 

dt 


= bpv cos zr 


(3) 


and 

A^^Ao+Ari+At-P (4) 

where Q = c-r e , A is the heat-transfer coefficient, 
p is the density of the free atmosphere, p m is the 
density of the meteoroids, c is the specific heat of 
the meteoroid, 6 is the reciprocal of the density 
scale height of the atmosphere, and the tempera- 
ture of the meteoroid outside the atmosphere is 
taken as zero. Equation (4) is simply an approxi- 
mate polynomial expansion of second degree in 
time for the heat-transfer coefficient. 

The following idealized history of each small 
fragment is assumed. The fragment leaves the 
parent body with temperature r c . Starting with 
this temperature, the fragment is heated to the 
evaporation temperature, cooling by radiation 
being taken into account. Isothermal heating is 
assumed for such small fragments, and the follow- 
ing equations are thus valid for the first part of 
the trajectory: 


C Z^^f+a R (T*-T ,«)- 


dv _ -3Tp,v 2 
dt 4 RfPm. 


II 

— A (t)pv* 
8 Qpm 

(1) 

dH 

= — V COS Zr 

(7) 

dv 

dt 

-3T P v 2 
4 Rpm 

(2) and 

dt 

P/ =P~F(H 0 -H) 

(8) 
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Here, <tr is the Stcfan-Boltzmann constant, A/ 
the heat-transfer coefficient for fragments, and 
T the drag coefficient. 

We deal with a very simplified case in this paper, 
the air density at height H being the only dis- 
turbed entity; the density in the wake is denoted 
by p f . In equation (8), F describes this dis- 
turbance. The function F also depends on the 
height of fragment separation, H 0 ; it has two 
extreme values : F = 0 if the atmosphere behind the 
parent body is not disturbed, and F~p if there 
is a vacuum behind. 

After being heated to the evaporation tem- 
perature, the fragment starts the second part of 
its trajectory. The temperature is determined by 
radiation, vaporization, and aerodynamic heating. 
Mass loss and luminosity occur until the energy 
flux is balanced by radiation. The second part of 
the fragment’s trajectory can be described by the 
following equation : 

dRj A /P/ ^4^(T 4 — 7V) 

dt 8 Q fPm ( } 


are not shining, because they are not yet heated 
enough. The beginning of the light curve is taken 
at the height corresponding to the air density p B 
computed for the ending of the preheating period 
according to the formula by Ceplecha and PadevSt 
(1961): 




8Xr(JZ) / e™+e-™ 

aV'JR \ e wR_ e -wR 



( 11 ) 


where W= (b cos zirv^) 112 / p } (3 2 is the thermal 
diffusivity, X is the thermal conductivity of the 
meteoroid, r is the temperature at which struc- 
tural failure occurs, and a is the accommodation 
coefficient. The end of the light curve is defined 
to be at the height where the velocity has de- 
creased to 5 km/s. 

The dynamic and photometric masses can now 
be computed and compared by means of the same 
procedures as for observed meteors. Equations 
(1) to (4) are used to find the radius and thus the 
mass of the body, and the conventional luminous 
equation 


and by equations (6), (7), and (8). Here, Q f is 
the heat of ablation of vapors from the fragment. 
Visual radiated luminosity If belonging to one 
fragment is then described in terms of the lumi- 
nous efficiency coefficient to by 


if=y 2 A(y S Tp m y!*Tov*R f \ Ho) 




X CfTpfV 


, 3 __ 


4 

Qf 


j 


X exp T- J (»//„*—»») (10) 


where 


A/ 

c ‘ 2r q , 

Also, A denotes the shape factor. 


PROCEDURE USED 

The total light intensity at a given height H m 
is the sum of all the partial intensities from all 
the fragments that separated from the parent 
body at heights H>H m . It is clear that the frag- 
ments that separated within heights Ih>H>H m 


2 f **** lit') 

m p {t) = - ~~ dt'+m end (12) 

To J t V 3 

is used to determine the photometric mass. Here, 
v is the velocity of the parent body and I is the 
instantaneous intensity. The integration constant 
Wend was taken as the dynamically determined 
mass at the terminal point. 

The change of air density behind the parent 
body might alter the computed photometric 
masses, but there is no need to attain higher 
values than those in the paper by McCrosky and 
Ceplecha (1970). However, if, for example, a 
constant factor for decreasing the density behind 
the main body is used, this undesired effect 
happens. But such a simple model for the function 
F docs not correspond even to a primitive physical 
guess. We would expect a decrease of the dis- 
turbance farther behind the meteor body, ter- 
minating with the density of the undisturbed 
atmosphere. It seems that each realistic model of 
the disturbance will result in a decrease of lumi- 
nosity originating from one fragment. In an 
attempt to increase the photometric mass, we 
should find a model of the disturbance that 
would decrease the level of the photometric mass 
as little as possible and we should look for another 
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Figure 1. — Logarithm of the ratio of the photometric 
mass nip to the dynamical mass titd for an initial 
radius = cm and an initial velocity tu = \2 
km/s. 



Figure 3. — Logarithm of the ratio of the photometric 
mass nip to the dy namical mass for an initial radius 

#*> = 100 cm and an initial velocity i’* = 30 km/s. 



Figure 2. — Logarithm of the ratio of the photometric 
mass nip to the dynamical mass mj for an initial radius 
Rn — 1000 cm and an initial velocity v * = 12 km/s. 
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effect strong enough to increase the level. The 
first condition may be met by a rather extreme 
and simple function — the discontinuous function 
F=p for an interval of H 0 >H>H 2 , and F = 0 for 
an interval of The fragments begin 

by moving in vacuum but cannot move ahead of 
the decelerating parent body. As an extreme case, 
we can imagine that the fragments are moving in 
a trap close behind the body, departing at height 
H 2 . This extreme choice of the function F serves 
the computations well as the first approximation 
of some more realistic model. The real function 
should start from a value very close to the air 
density (almost vacuum), followed by a steep 
change that is represented here by the discon- 
tinuous jump to zero. All our preliminary com- 
putations used this discontinuous function. 

An additional parameter now is the number of 
fragments crossing the discontinuity in the air 
density at a given height H m . If the discontinuous 
function F is as above, only one change results: 
The height of the separation of the corresponding 
fragment will be higher than in the case without 
any disturbance. But we are interested in differ- 
ences of orders of magnitude, and this is not 
possible with the above model with A constant. 
On the other hand, if we assume that A decreases 
with time, a steep decrease of the main-body mass 
results during the second part of the trajectory. 
Generali}^ speaking, the decrease of A with time 
can be expected from the physical point of view". 

RESULTS 

Until now', only a limited number of cases have 
been computed. The numerical choice of param- 
eters w r as made rather arbitrarily, with only great 
differences between the computed photometric 
and the computed dynamic masses being ac- 


cepted. We can guess that the differences will be 
less after more physical factors are includod, but 
the computed examples illustrate the general 
situation well. 

Using the above model, w r e can explain the 
great difference between the dynamic and the 
photometric mass during the second part of the 
trajectory, if w r e assume chondritie densities for 
the body. As a function of time, the difference 
between the two masses is less in the first part of 
the trajectory and becomes greater at the end 
(see figs. 1 to 4) . 

COMMENTS 

Z) = Bq-{- H 2 -\- B 2 • Z/ 2 2 is the difference in 
height in kilometers between the fragment 
separation (H 0 ) and that of departure from the 
trap (H 2 ), 

Numerical parameters used for computation 
of the values in tables 1 to 4 (cgs) are as follows: 


a 

= 1 

A 

= 1.21 

b 

= 1.6 X 10~ 6 cm -1 

c 

= 10 7 cm 2 s -2 deg -1 

Qr 

= 8X 10 10 cm 2 s~ 2 

initial Rj 

r = 0.01 cm 

To 

w 

o 

o 

II 

T b 

= 2580° K 

COS Zr 

= 1 

r 

= 0.5 

A 

= 3 X 10 5 g cm s -3 deg -1 

Pm 

= 3.5 g cirr 3 


A/ 

= tv, — 2X10 -12 cm _2 s 2 

T c 

= 600° C 


The results do not depend on the value r 0 . 
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Velocity v 
(km/s) 


5.000 
5.692 
6.446 
7.225 
7.980 
8.669 
9.268 
9.772 

10.191 

10.536 

10.822 

11.058 

11.255 

11.421 

11.561 

11.680 

11.780 

11.866 

11.939 

12.000 


mamical mass mj 
(g) 

Height H 
(km) 

Epoch t 
<■) 

Assumed 
heat-transfer 
coefficient A 

6.0967X10* 

17.654 

4.2592 

0.00224 

6.4917X10* 

17.999 

4.1959 

0.01734 

7.8346X 10 2 

18.363 

4.1370 

0.03138 

1 . 0994 X 10 3 

18.750 

4.0816 

0.04458 

1 .7900X10* 

19.161 

4.0286 

0.05720 

3.2799 X10 3 

19.601 

3.9768 

0.06952 

6.4496X10 3 

20.072 

3.9253 

0.08178 

1.3035X10 4 

20.570 

3.8728 

0.09426 

2.6252X10 4 

21.111 

3.8184 

0.10719 

5.1745X10 4 

21.706 

3.7610 

0.12082 

9.9022X10 4 

22.365 

3.6994 

0.13545 

1.8352X10 5 

23.104 

3.6322 

0.15140 

3.2981 X10 5 

23.945 

3.5576 

0.16911 

5.7708X 10 5 

24.920 

3.4727 

0.18922 

9.8699X10 5 

26.079 

3.3737 

0.21267 

1 . 6587 X 1 0 6 

27.506 

3.2540 

0.24100 

2.7585X10* 

29.358 

3.1011 

0,27714 

4.5839X10 6 

31.994 

2.8875 

0,32754 

7.7595X10* 

36.419 

2.5258 

0,41270 

1.4661 X10 7 

73.536 

0.0000 

1,00000 


Velocity v 
(km/s) 


5.000 
6.281 
7.454 
8.408 
9.145 
9.708 

10.144 

10.487 

10.762 

10.987 

11.174 

11.330 

11.463 

11.576 

11.674 

11.759 

11.832 

11.896 

11.952 

12.000 


Table 2. — Meteoroid With Initial Velocity v^-12 km/s and Radius R * - WOO cm 


Dynamical mass 
(g) 


Height ff 
(km) 


4 . 3224 X 10 4 
8.8988X1G 4 
2.4031 X10 5 
7.1522X10 5 


1.7160 

2.2542 

2.8161 

3.4026 


2.0631 XlO 6 
5 . 4479X10 6 
1.3075X10 7 
2.8694X10 7 
5.8364X10 7 
1.1141X10 8 
2.0175X10* 
3. 4984X10* 
5.8537X10* 
9.5153X10* 

1 .5114X10® 
2. 3604X10® 
3.6457X10® 
5.6144X10® 
8.7454X10® 
1.4661 X10 10 


4.0220 
4 , 6686 
5.3541 
6.0834 
6.8585 
7.6891 
8.5869 
9.5653 
10.643 
11.699 
12.860 
14.281 
16.114 
18.698 
23.053 
73.532 


Epoch t 

(a) 


Assumed 

heat-transfer 


coefficient A 


5.2882 
5.2282 
5.1762 
5.1285 
5.0826 
5 . 0370 
4.9903 
4.9419 
4.8907 
4.8360 
4.7766 


0.05170 

0.06718 

0.08050 

0.09263 

0.10425 

0.11574 

0.12742 

0.13949 

0.15215 

0.16561 

0.18011 


4.7112 0,19596 

4.6380 0.21355 


4.5544 

4.4564 

4.3374 

4.1851 

3.9716 

3.6087 

0.0000 


0.23343 

0.25647 

0.28406 

0.31878 

0.36624 

0.44380 

1.00000 
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Table 3. — Meteoroid With Initial Velocity v w =S0 km/s and Radius 'R ao =100 < 


Velocity v 
(km/s) 


5.000 

23.680 

26.363 

27.493 

28.144 

28.571 

28.874 

29.101 

29.276 

29.416 

29.530 

29.623 

29.701 

29.766 

29.822 

29.870 

29.910 

29.945 

29.975 

30.000 


Dynamical mass ms 
(g) 

Height H 
(km) 

Epoch t 
(«) 

3.2649X10 -6 

29.104 

1.7436 

7.7116X10-' 

29.455 

1.7215 

2.8707X10' 

29.826 

1.7074 

2.2584X10 3 

30.220 

1.6934 

9.9065X10" 

30.641 

1.6789 

3.1699X10 3 

31.091 

1.6637 

8.3382X 10 3 

31.576 

1.6476 

1.9233X10 4 

32.096 

1.6303 

4.0281 X10 4 

32.646 

1.6117 

7.8608X10 4 

33.250 

1.5914 

1.4512X10 5 

33.919 

1.5691 

2.5621 X 10 5 

34.675 

1.5442 

4. 3758 X 10 s 

35.540 

1.5161 

7.2682X10 5 

36.548 

1.4837 

1 . 1836 X 10 s 

37.757 

1.4455 

1 .9020X10® 

39.260 

1.3988 

3.0356X10 6 

41.235 

1.3387 

4.8632X10 6 

44.094 

1.2542 

7.9709X 10 6 

49.320 

1.1101 

1.4661 X10 T 

87.079 

0.0000 


Assumed 
heat-transfer 
coefficient A 


0.01609 

0.02887 

0.03705 

0.04514 

0.05352 

0.06230 

0.07159 

0.08154 

0.09226 

0.10396 

0.11683 

0.13114 

0.14733 

0.16591 

0.18784 

0.21462 

0.24899 

0.29726 

0.37931 

1.00000 


Table 4. Meteoroid With Initial Velocity v ,=30 km/s and Radius R„ = 1000 cm 


Velicity v 
(km/s) 

Dynamical mass 
(g) 

Height H 
(km) 

Epoch t 
(s) 

Assumed 
heat-transfer 
coefficient A 

5.000 

8. 6440X10-’ 

13.510 

2.2721 

0.01375 

24.240 

2.6803X10 3 

13.855 

2.2333 

0.03189 

26.458 

6.9058X10 4 

14.219 

2.2193 

0.03842 

27.473 

4.5678X10 5 

14.605 

2.2052 

0.04498 

28.083 

1 .7743X10® 

15.016 

2.1907 

0.05173 

28.496 

5.1904X10 6 

15.455 

2.1755 

0.05881 

28.799 

1.2789X10 7 

15.928 

2.1593 

0.06634 

29.030 

2. 7879X10 3 

16.438 

2.1420 

0.07438 

29.211 

5.5661 X10 7 

16.993 

2.1234 

0.08303 

29.358 

1 .0420X10 8 

17.600 

2.1030 

0.09246 

29.479 

1 .8520X10* 

18.272 

2.0807 

0.10280 

29.579 

3. 1676 X 10 s 

19.024 

2.0557 

0.11433 

29.664 

5.2456X10 8 

19.875 

2.0276 

0.12732 

29.736 

8. 4820 X 10 s 

20.838 

1.9951 

0.14225 

29.798 

1. 3466X10 9 

21.976 

1.9568 

0.15984 

29.851 

2. 1108X10 8 

23.377 

1.9101 

0.18126 

29.897 

3.2897X10® 

25.198 

1.8499 

0.20871 

29.937 

5.1470X10® 

27.790 

1.7652 

0.24718 

29.971 

8.2266X10® 

32.286 

1.6207 

0.31234 

30.000 

1 .4661 X10 10 

87.077 

0.0000 

1.00000 



